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Does a single eigenstate encode the full Hamiltonian? 
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The Eigenstate Thermalization Hypothesis (ETH) posits that the reduced density matrix for a 
subsystem corresponding to an excited eigenstate is “thermal.” Here we expound on this hypothesis 
by asking: for which class of operators, local or non-local, is ETH satisfied? We show that this 
question is directly related to a seemingly unrelated question: is the Hamiltonian of a system encoded 
within a single eigenstate? We formulate a strong form of ETH where in the thermodynamic limit, 
the reduced density matrix of a subsystem corresponding to a pure, finite energy density eigenstate 
asymptotically becomes equal to the thermal reduced density matrix, as long as the subsystem size 
is much less than the total system size, irrespective of how large the subsystem is compared to 
any intrinsic length scale of the system. This allows one to access the properties of the underlying 
Hamiltonian at arbitrary energy densities/temperatures using just a single eigenstate. We provide 
support for our conjecture by performing an exact diagonalization study of a non-integrable ID 
lattice quantum model with only energy conservation. In addition, we examine the case in which 
the subsystem size is a finite fraction of the total system size, and find that even in this case, a 
large class of operators continue to match their canonical expectation values. Specifically, the von 
Neumann entanglement entropy equals the thermal entropy as long as the subsystem is less than half 
the total system. We also study, both analytically and numerically, a particle number conserving 
model at infinite temperature which substantiates our conjectures. 
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I. INTRODUCTION 

Given a local Hamiltonian, what information about the 
system is encoded in a single eigenstate? If the eigenstate 
happens to be a ground state of the Hamiltonian, tremen¬ 
dous amount of progress can be made on this question for 
Lorentz invariant system^ 1 2 , especially conformal field 
theories (CFTs)P®, and for topological phased®^. For 
example, one can read off the central charge of a CFT 
from the ground state entanglement®^, while for topo¬ 
logical phases, essentially all Topological data’ such as 
braiding statistics of anyons can be extracted from the 
degenerate ground states 8 10 . In this paper we argue 
that a single finite energy density eigenstate of an er- 
godic quantum many-body Hamiltonian is sufficient to 
determine the properties of the system at all tempera¬ 
tures. 

It is not very surprising that the ground states of quan¬ 
tum many-body systems contain some information about 
their excitations. This is because an entanglement cut 
often mimics an actual physical cut through the system, 
thus exposing the underlying excitations along the entan¬ 
gling boundary®. The same intuition is tied to the fact 
that the ground state entanglement satisfies a “boundary 
law” of entanglement entropy 11 ^, that is, the von Neu¬ 
mann entanglement entropy Si = — tr^ (pa log(pA)) of 
the ground state corresponding to a subsystem A scales 




2 


with the size of the boundary of subsystem A. 

How does the nature of information encoded evolve as 
one goes from the ground state to an excited eigenstate? 
Typically, there always exist eigenstates with energy E 
just above the ground state which continue to satisfy 
an area law of entanglement. These are the eigenstates 
which have zero energy density, i.e. limy^oo E ^fi° = 0 
where Eq is the ground state energy and V is the total 
volume of the system. These eigenstates can often be in¬ 
terpreted as the action of a sum of local operators acting 
on the ground state; for example, in a system with spon¬ 
taneous symmetry breaking one can construct an eigen¬ 
state consisting of a few magnons by a superposition of 
spin-flips acting on the ground state. Furthermore, the 
level spacing between two contiguous low-lying excita¬ 
tions scales as SE ~ 1 /L a where a > 0 depends on 
dimensionality and the phase of matter under consider¬ 
ation. In this paper, we will instead be concerned with 
excited eigenstates that have a finite energy density , i.e. 
limy^oo E ^fi° fiz 0. For notational convenience, we will 
set Eq = 0 for the remainder of this paper. 

As argued by Srednickf 13 [ a typical finite energy den¬ 
sity state (i.e. a typical state in the Hilbert space that 
satisfies = Ve where e is the energy density) 

when time-evolved with the Hamiltonian H for suffi¬ 
cient time is expected to lead to predictions dictated by 
the basic tenets of equilibrium statistical mechanics, if 
the system thermalizes. Such an expectation leads to 
the “Eigenstate Thermalization Hypothesis” (eth)EMH] 
which stipulates that the thermalization occurs at the 
level of each individual eigenstate. An alternative ap¬ 
proach by DeutsclP^, which is based on perturbing an 
integrable system by a small integrability breaking term, 
leads to the same suggestion. If ETH holds true, then in 
the thermodynamic limit the equal-time correlators of an 
operator with respect to a finite energy density eigenstate 
\fi) are precisely equal to those derived from a thermal 
ensemble, i.e. 




tr (Oe f3H ) 
tr ( e~P H ) 


(1) 


where /3 is chosen such that the Eq. [l] holds true when 
O = .H, the Hamiltonian. Henceforth we will use the no¬ 
tation | to denote an eigenstate whose energy density 
corresponds to temperature /3 -1 . A notable exception 
to ETH is a many-body localized system in the context 
of strongly disordered interacting quantum systems, 16 22 
which fails to thermalize and does not satisfy Eq. [l] The 
possibilit)BMI3 or impossibility®® of the violation of 
ETH without disorder has also been discussed recently. 

In this paper, we restrict ourselves to systems where 
ETH, as defined by Eq. [l] holds. However, Eq. [l] alone is 
incomplete unless one also specifies the class of operators 
for which it holds. For example, one simple non-local 
operator for which Eq. [l] breaks down is the projection 
operator onto the eigenstate \fi) that enters Eq. 

[lj the left hand side of Eq. [l] yields unity for this oper¬ 


ator, while the right hand side is exponentially small in 
the volume, a clear disagreement. On that note, it is of¬ 
ten mentioned that in systems where Eq. [I] does hold, it 
does so only for “few body” op er at or P 4 36 Iwhere, to our 
knowledge, the precise meaning of few-body operator has 
not been clarified (see Ref. [37] for related discussion). In 
this paper, we conjecture and provide numerical evidence 
that Eq. [l] holds for all operators within a subsystem A 
when the volume Va of subsystem A satisfies < F 
(or, more precisely, when Va/V — 0 as V oo). We 
also explore the more general case where subsystem A 
spans a finite fraction / = Va/V > 0 of the total system 
size. We provide some evidence that when the fraction is 
less than a critical 0(1) number /*, then all operators not 
explicitly involving energy conservation take their ther¬ 
mal values. We also explore the more general condition 
Va < V/2 and show that even in this case, Eq. [l] holds 
for a large class of operators. On that note, we should 
mention that the questions such as which Hamiltonians 
(and which operators) satisfy ETH is now entering the 
realm of experimental physics (see e.g. Ref. l38l) due to 
advances in high resolution imaging techniques^. 

The satisfaction of Eq. [l] for all operators in a subsys¬ 
tem A is equivalent to the statement that the reduced 
density matrix pa(\^)p) = tr^^)^^! corresponding to 
an eigenstate \fi)p is given by 


Pa{ H)p) = PaMP) O) 


where 


PA, th(/3) 


fa*(e-* g ) 

tr (e~P H ) ’ 


A being the complement of A. Note that the trace in 
the denominator is over the whole Hilbert space. When 
Va is held constant, the equality in Eq. [2a] means the 
density matrices become elementwise equal in any basis 
as V oo. 

One immediate consequence of Eq. [2a] is that the ther¬ 
modynamical properties of a system at arbitrary temper¬ 
atures can be calculated using a single eigenstate. For 
example, Eq. [2a] implies that to the leading order, the 
Renyi entropies S a (= — log [tr^(p^)]) for an eigen¬ 
state \fi)y corresponding to a subsystem A with Va <C V 
are given by 


S a = -^—AaP (/(<*/?) - /(/?)) , (3) 

a — 1 

where /(/?) is the free energy density at temperature /3 _1 . 
The above equation allows one to access the free energy 
density / at an arbitrary temperature by varying a. Note 
that Eq. [3] holds only to the leading order because Renyi 
entropies S a receive additional subleading contributions 
due to the conical singularity induced at the boundary 
of subsystem A 3 5 . In the limit a 1, one recovers 
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the equality between the von Neumann entanglement en¬ 
tropy Si and the thermal entropy 5 t h = VASth(P), where 
s th{/3) is the thermal entropy density at temperature /3 _1 , 
a result which was argued to hold in Ref. TO for the spe¬ 
cial case of two weakly coupled ergodic systems. We em¬ 
phasize that these results cannot be derived from Eq. [l] 
alone were it to hold only for local operators, since entan¬ 
glement entropies do not correspond to the expectation 
value of any local operator. We also note that Refs. 1411421 
simulated the thermal Renyi entropy S a (starting with 
the expression on the right hand side of Eq. |2a| ) using 
Quantum Monte Carlo to access the properties of the sys¬ 
tem at temperature ( a/3 )~ 1 . Of course, Quantum Monte 
Carlo methods are not well suited to verifying ETH since 
they cannot access properties of a single eigenstate (the 
left hand side of Eq. 2a). 


We will also discuss an approximate, but more intuitive 
form of ETH, given by 


Pa( W)p) 


e -BH A 


tr a (e~P HA ) 


(2b) 


where Ha is the projection of the original Hamiltonian 
onto subsystem A. This form is approximate compared 
to Eq. [2a] because generically, it does not capture the cor¬ 
relations near the boundary correctly due to the some¬ 
what arbitrary truncation scheme used to obtain Ha • 
Nevertheless, equations [2a] and [2b] both yield the same 
results for all bulk quantities such as the Renyi entropy 
densities, as well as correlation functions of operators 
that have support only far from the boundary. 

A central task of this paper is to check the validity of 


Eqs. [2a] and [2b] and their consequences for model non- 
integrable systems. As already mentioned, we will argue 
that ETH allows one to calculate thermodynamical quan¬ 
tities as well as correlators at all temperatures/energy 
densities using only a single eigenstate. We will demon¬ 
strate this explicitly by studying a quantum ID model 
numerically. 

As mentioned above, we find evidence that Eq. |T] holds 
for many operators even when Va/V is held constant 
with Va/V less than some number /* > 0. In particu¬ 
lar, as we discuss later, our results strongly indicate that 
/ < 1/2 is sufficient to guarantee equivalence between the 
von Neumann entropy density of a pure eigenstate, and 
the thermal entropy density at the corresponding temper¬ 
ature. This is in contrast to Ref. ssi where it was argued 
that such an equivalence holds only in the limit /* —>■ 0. 
Recent 1}@E^, the requirement /* —>> 0 was substantiated 
using analytical and large scale numerical calculations for 
free fermions, an integrable system. Our results indicate 
that the /* —>> 0 requirement is likely a consequence of 
the integrable nature of the models in Refs. 1441451 

The paper is organized as follows. Sec. |TT| discusses 
general considerations for the validity of ETH, and in¬ 
troduces a division of all operators in a given subsystem 
into two distinct classes, which have different require¬ 
ments for ETH to hold. Sec. 2 ! illustrates some general 
features of ETH by studying properties of a hardcore bo¬ 


son model with global particle number conservation for 
infinite temperature eigenstates. Sec. [IV] introduces the 
model we study in the remainder of the paper, the trans¬ 
verse field Ising model with longitudinal field. Sec. [V] 
focuses on the entanglement entropies at finite tempera¬ 
ture. Sec. |VI| studies the validity of ETH when Va <C V 
by providing a close look into the entanglement Hamilto¬ 
nian, focusing on its spectrum and Schmidt vectors. This 
section also demonstrates the validity of Eq. [2a] when 
Va V by considering the trace norm distance between 
both sides. Sec. |VII| explores the validity of ET H when 
Va/V is taken to be finite as k —)> oo. Sec. VIII| provides 
an application, by using the reduced density matrix from 
a single eigenstate to predict correlators at all (finite) 
temperatures. Sec. m summarizes our results and pro¬ 
vides thoughts for future discussion. 


II. GENERAL CONSIDERATIONS 

A. Determining the Hamiltonian from microstates 
in classical statistical mechanics 

Suppose, for an isolated system described by classical 
statistical mechanics in a total volume V, we are given 
access to all classical microstates in a small energy win¬ 
dow [E,E + AE\, where A E ~ y/V is on the order of the 
energy fluctuations in the total system were the system 
coupled to a thermal bath, and thus all microstates corre¬ 
spond to the same energy density. We pose the question: 
does this information suffice to determine the underly¬ 
ing Hamiltonian, assuming that the Hamiltonian is local? 
The answer is indeed yes, following the standard proce¬ 
dure of obtaining the canonical ensemble from a micro- 
canonical ensemble. In particular, let us make a fictitious 
division of the system into A and A such that Va V^, 
and count the number of times a particular configura¬ 
tion Ca appears in subsystem A. This determines the 
probability distribution for finding a given configuration, 
P(Ca )• If all microstates are equally likely, then® 

E(C a ) 

P{Ca) = H { c A} e-mcA) (4) 

where E{Ca) is the energy in subsystem A. One may 
now invert this equation to obtain the energy E(Ca) = 
— ^ log(P(CT)), up to an irrelevant constant shift of en¬ 
ergy. In a classical statistical mechanical system E(Ca) 
is the Hamiltonian for subsystem A. In particular, know¬ 
ing E(Ca ), one may now calculate any thermodynamic 
property at any temperature. Here it is crucial to note 
that Eq. [4] does not assume that the energy density 
E(Ca)/Va equals the energy density E/V of the mi¬ 
crostates being sampled. 

As discussed in the introduction, we will provide ev¬ 
idence that the quantum mechanical analog of Eq. [4] is 
given by Eqs. [2a] |2b| We now proceed to discuss the 
conditions under which Eqs. [2a] [2b] are valid. 
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B. Two classes of operators 


For reasons soon to be discussed, we find it useful 
to separate operators in a given Hilbert space into two 
classes: 

Class I (“Equithermal Operators”): If the reduced 
density matrix takes the thermal form (i.e. the right hand 
side of Eq. 2b), then in the limit Va oo, the expecta¬ 


tion value of equithermal operators receives contribution 
only from the eigenstates of H A at an energy density 
corresponding to the temperature /3 _1 . One might have 
thought that this is true for all operators, however, there 
exist operators such as e ~ n ^ HA , whose expectation value 
includes contribution from eigenstates of Ha at temper¬ 
ature [(n + 1)/3] _1 in addition to the temperature f3~ 1 . 
Clearly, local operators fall into this class, as do sums 
of local operators. Several non-local operators, including 
the von Neumann entropy Si, also fall into this class. 

Class II ( u Non-equithermal Operators”): We dub 
all operators not in Class I as “non-equithermal opera¬ 
tors”, or Class II operators. All Renyi entropies S a (for 
a^l) fall into this clasiP^. 


C. ETH: Class I vs. Class II operators 

Let us first consider the relationship between Eq.[l]and 
Eqs. [2a| |2b| Eq. 1 may be rewritten as, 


tr A ( PaO ) = 


tryt(Qtr^(e P H )) 
tr ( e~P H ) 


(5) 


If this equation holds for all operators in a subsystem 
A, hermitian as well as non-hermitian, then one obtains 
Eq. ^a| pa(\^)p) = Pa, th(/3)- This is because one may 
expand both the p A and p A ,th i n terms of the complete 
set of operators in subsystem A , and by choosing appro¬ 
priate O prove that they are equal to each other element- 
by-element. One of the most important consequences of 
this equality is that it allows one to extract properties of 
the Hamiltonian at arbitrary temperatures using a sin¬ 
gle eigenstate, which is one of the central points of this 
paper. 

We will now discuss ETH for both Class I and Class 
II operators. For each class, we consider separately two 
cases: (i) when Ih < L; and (ii) when the ratio / = 
Va/V is taken to be fixed and finite as Va, V oo. 


1. ETH for Class I operators 

Let us first briefly discuss ETH for Class I operators 
when Va <C V. This includes both the case where V A is 
held fixed as F oo, as well as the case where the limits 
Va,V oo are taken such that V A /V 0. In fact, 
this is the traditional definition of ETH—that all local, 
“few-body” operators match their values in the canonical 
ensemble in this case. 


Let us now consider the validity of ETH for Class I 
operators in the fixed-ratio limit where 0 < / < | is 
finite. In contrast to classical statistical mechanics, we 
expect that quantum mechanically, one does not require 
the constraint V A Vj for ETH to hold for a large class 
of Class I operators. Indeed, as discussed below, several 
known results already point to the conclusion that Eq. [l] 
holds for at least some operators in Class I, as long as 
Va < Vj with both Va, Vj 00 • 


One piece of evidence that suggests that Eq. [l] might 
hold for Class I operators as long as Va < Vj comes from 
the study of quantum quenches in conformal field theo¬ 
ries (CFTs). As shown in Ref.[48j the time-dependent re¬ 
duced density matrix p^(t) of a system initially prepared 
in a low-entanglement state, and evolved with a CFT 
Hamiltonian, approaches the thermal density matrix, as 
long as Va < E/2, with Va, V oo. Ref. [H] character¬ 
ized the closeness between p A {t) and the thermal density 
matrix p A ,th (Eq. 2a) in terms of the operator overlap 
lit) = t— tr (pa (t)p A ,th) ^ which is exponentially close 

(tr(p 2 A (t))tr(p 2 At J) 

to unity for Va/2 < t < Vj/2. It is important to note 
that in the thermodynamic limit, I only receives contri¬ 
bution from eigenstates at temperature /3 _1 , so this only 
guarantees that operators in Class I will satisfy Eq. [l] 

Another piece of evidence comes from the recent stud¬ 
ies of large central charge conformal field theories^ 51 . 
In particular, Refs. I49|51l studied the entanglement en¬ 
tropy of pure eigenstates in finite temperature confor¬ 
mal field theories with large central charge. In the limit 
Va, V 1/T, while keeping V A /V fixed, it was found 
that the entanglement entropy becomes equal to the 
thermal entropy at all non-zero temperatures as long as 
Va < Vj. 

Lastly, the entanglement entropy for a random pure 
state is given 


5=-iog(i^r i + i^r i -i^r 1 ) (6) 

where \Ha\, 1^1 are s i zes °f the Hilbert spaces 

of subsystems A , A and the total system (= A U A) 
respectively. Thus, as soon as Va < Vj, one obtains 
S = — logd'H^I), which is indeed the thermal entropy 
for subsystem A at infinite temperature. Since random 
pure states mimic eigenstates at infinite temperature (i.e. 
|V ; )/3 =o) 5 this again suggests that the condition Va < Vj 
is perhaps sufficient, at least for some operators. 

On the other hand, there is a well-known Class I op¬ 
erator for which ETH fails when the ratio of subsystem 
to total system size f = V A /V is finitd^. When / is 
finite, the energy variance of the reduced density ma¬ 
trix Pa(\^)b) suppressed by a factor of (1 — /) 

compared with the value the variance would have taken 
in the canonical ensemble. Ultimately, this is due to the 
fact that a single eigenstate has precisely zero energy vari¬ 
ance ((H — (H)) 2 ) in the full system, unlike the canonical 
ensemble, where the variance scales proportionally with 
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system size. This relationship can be expressed as 

tr[p A {Wp)0 A ,p\ = jr te[pA,th(P)°A,p], (7) 

where 0 A ,p = {Ha — {H A ) p) 2 is the energy variance op¬ 
erator. We will explore implications of the subsystem 
energy variance mismatch more carefully in Sec. m 

It is worth noting that for the case of time-evolved 
states, the full system variance is independent of time. 
For a given initial state, this variance may indeed be dif¬ 
ferent from the energy variance expected in the canonical 
ensemble, which implies that the energy variance for any 
subsystem that is a finite fraction of the total system will 
disagree, even at long times 56 . However, we expect that 
for “typical” initial states (which are typically inaccessible 
from an experimental point of view), the overall energy 
variance will match its result in the canonical ensemble, 
and so the energy variance for any subsystem will also 
match after thermalization (“canonical typicality’®^) . 

Overall, while we expect that ETH is obeyed by many 
Class I operators when / > 0 is finite, it cannot be 
satisfied by all such operators, since the subsystem en¬ 
ergy variance provides an important counterexample. 
Nonetheless, we expect that all Class I operators not re¬ 
lated to energy conservation (or another conserved quan¬ 
tity) should satisfy ETH as given in Eq.[l] A more precise 
conjecture along these lines is that the set of operators 
spanned by H\ where n ranges between unity and the 
size of the Hilbert space do not satisfy ETH in the sense 
of Eq. [7] above. The number of such operators is still ex¬ 
ponentially smaller than the total number of independent 
operators in a subsystem (e.g. in a spin-1/2 spin system, 
the total number of operators in a region A is 4 Va , while 
the number of operators of the form H\ is 2 Va , the size 
of the Hilbert space in region A). 


2. ETH for Class II operators 


ETH will hold for all Class II operators either when / 
is finite. In addition, for a given ratio Va/V with both 
Va,V —>• oo, there is a physical constraint on the range of 
energy densities for which the spectrum of \'if)p in prin¬ 
ciple can match that of p A ,th(f3 )• To appreciate this, 
let us consider a slightly different problem—an arbitrary 
Hamiltonian of hardcore bosons with particle number 
conservation, at infinite temperature. We will consider 
an explicit example of such a system in the next sec¬ 
tion. Since the total particle number operator N com¬ 
mutes with the Hamiltonian and satisfies the equation 
N = Na + N A7 the reduced density matrix pA for a 
wavefunction |- 0)^=0 is block diagonal in the number of 
particles Na in subsystem A. Furthermore, if ETH holds 
(as given by a generalization of Eqs. 

Schmidt decomposition is given by 

N 

IV’)/3=0 = E AE E I u i)N A ® \Vj)N-N A ( 8 ) 

N a =0 i 

where \n a are the Schmidt coefficients in the sector Na, 
and \ui) n a , \vf) n-n a are the corresponding eigenvectors. 
The label i captures fluctuations of particles within a 
fixed sector Na- Note that there is no index i on Xn a 
because we are at infinite temperature and all Schmidt 
states within a sector Na are equally likely. 

The decomposition in Eq. [5] allows one to calculate 
properties of subsystem A at infinite temperature even 
away from filling N/V since the reduced density matrix 
Pa will contain sectors with various densities N A /V A - 
However, there is both an upper limit and a lower limit 
on the density in subsystem A , since 


2a and 2b), then the 


max [N — {V — V A ), 0] < N a < min [AT, V A ] (9) 

And thus the particle density Na/Va in subsystem A 
satisfies 


The extra ingredient introduced by Class II opera¬ 
tors is that if ETH holds for them, then taking such 
an operator’s expectation value with respect to a state 
IV0/3 allows one to access the properties of the Hamilto¬ 
nian at a temperature different than /3 _1 . For example, 
the Renyi entropy S a corresponding to p A {\^)p) satis¬ 
fies S a = p^Va/ 3 (f(af3) — /(/?)), thus allowing one to 
access the free energy density at temperature {a/3)~ 1 . 

Let us first consider the validity of ETH for Class II op¬ 
erators when V/ Ck. Remarkably, the results presented 
in the remainder of this paper demonstrate that ETH is 
valid for all Class II operators in this limit. Thus, a single 
eigenstate of finite energy density contains knowledge of 
the properties of the system at all temperatures. 

Now let us turn to the case in which Va/V is finite, 
which turns out to be much more subtle. As mentioned 
in the previous subsection, there is already a Class I op¬ 
erator for which ETH fails in this limit, namely the sub¬ 
system energy variance. Thus, we do not expect that 


max [1 — (1 — n)//, 0] < ——- < min [n//, 1] (10) 

Va 

where n = N/V is the overall particle density and / = 
Va/V. Thus, a necessary condition for the wavefunction 
in Eq. [8] to encode properties of the system at all fillings 
is 

/ < min [n, 1 — n\ (11) 

The above discussion, with some modifications, carries 
to systems with (only) energy conservation, at an arbi¬ 
trary temperature. The Schmidt decomposition of an 
eigenstate 1 if)p with eigenvalue E may now be written 
as: 

IV>)/3 = J2 Al\ u i) ® k> ( 12 ) 
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The physical content of ETH, as approximated in Eq. 
2 b is that A^ oc e ~ (3EA ’ i where Ea,i is the V th energy 


eigenvalue of Ha (the projection of the Hamiltonian to 
subsystem A) while |iq) is the corresponding eigenstate 
of Ha • Denoting the ground state energy to be zero, 
one naively expects that (ui\HA\ui) < E V |iq) since the 
energy density in the subsystem A cannot be less than 
the ground state energy density. However, this argument 
has a loophole since in contrast to the particle number 
operator TV, the total Hamiltonian is not separable into 
subsystems A and A: H = Ha + H -4 + H AA , which 
actually allows (ui\HA\ui) to exceed E as we will see in 
Sec. |VII| in the context of the model Hamiltonian in Eq. 


30 below. To understand the constraint on (ui\HA\ui) 


precisely, let us derive an expression which encapsulates 
the classical notion that the sum of energies in subsystem 
A and A equals E. 

We first note: 


{uio\® (v i 0 \H\ip)p = E{u i0 \ <g> (vio\ip)p (13) 

= E^/X~ 0 (14) 

The above expression can be re-evaluated using the 
decomposition H = Ha + H -^ + H a a'- 

(«io| 0 (v i0 \H\ip)p (15) 

= (u i0 \0 (v i0 \H a + H a + H aa \ ip)p 

(‘M’iO | A | U i0 ) “l - \J A (0 \ | H~ A | ViQ ) 

E A'^ol ® {v i 0 \H AA \uj) <g> \vj) (16) 

3 


Equating the two ways to calculate the same expres¬ 
sion, one finds: 


{v i0 \H A \v i0 ) + E V ® ( Vio\H AA \uj) 0 | vj) 


E (uio | Ha |^io) 


(17) 


Due to the variational principle for the ground state, 
(vio\H A \vio) > —cL d ~ x where c is a constant (recall 
that in our convention, the ground state energy for the 
full Hamiltonian is set to zero). Since both E and 
(uio \ Ha |^io) scale as L d , the only way for (iqo|JC 4 |^o) to 
exc eed E is that the second term on the left hand side of 


Eq. 


17 


viz. -^boundary V . \j (^i 0 | ® (^iO\H- aa\^ j ) * 


| Vj), is negative and scales as L d . When that hap¬ 
pens, ETH no longer holds, as we now argue on general 
grounds, and will also demonstrate numerically for a lat¬ 
tice Hamiltonian in Sec. |vn] To see this, we reiterate that 
ETH requires that (i) |iq)’s are approximate eigenstates 
of Ha , and (ii) A^ oc e ~P( u i\ H A\ui) _ e -pE A ,i m Firstly, 
when (uio\HA\uio) < E so that ETH could in principle 
hold, the ^boundary term can be neglected because the 
‘diagonal term’ in ^boundary (ho. the term corresponding 
to j = iO) scales as the boundary (oc L d_1 ) and is thus 


subleading, while the off diagonal terms scale as e~ L 
and thus vanish in the thermodynamic limit (recall that 
V A > Va). °ther hand, when (i^o^aI^o) > E, 

the l^o)’s now correspond to states of zero energy den¬ 
sity, and the aforementioned argument for neglecting off- 
diagonal terms is no longer valid. So, let us assume that 
(uio\HA\uio) > E and each |iq 0 ) continues to be an eigen¬ 
state of Ha • Thus, one requires that 



A(eQ 

A(e) 


M(e, e')e s(e,) 


oc g(e)/L d 1 , 


(18) 


where we have taken the continuum limit and A(e) 
denotes the Schmidt eigenvalue corresponding to an 
eigenvector | u) at energy density e, while M(e, e') = 
(u(e)\ 0 {v(e)\H AA \u(e f )) 0 \v(e')) and g(e) = e - 


{u{e)\HA\u{e))/L d . It is obvious from Eq. 18 that A(e) oc 


e -/3E A _ e -/3efL d - 1S nQ } on g er the solution. In fact, the 
only way for the integral on the left hand side of Eq. [18] 
not to have any exponential dependence on L (as required 
by the right hand side) is that the integrand itself does 

not have such dependence, i.e. yj oc e /) e~ s( ^ e '\ 

This implies a breakdown of ETH when (ui$\H aIu^) > 
E. 


The above discussion implies that for a given wave- 
function and bipartition, the maximum energy density 
that is potentially accessible in a subsystem A , such that 
the corresponding Schmidt weight satisfies ETH is, 


e* = min (E/V A , e max ) = min(e//, e max ) (19) 

where e = E/V is the energy density corresponding 
to the wavefunction and e max is the maximum energy 
density for the Hamiltonian H (recall that e max can be 
finite for lattice-regularized quantum systems, e.g. for 
models of fermions or spins/hardcore bosons). Above, 
we have assumed that e < e max /2. In the case when 
e > e max /2, the range of available energies is instead 
bounded from below by max [0, e max (l — 1//) — e/f]. If 
our goal is to capture the fluctuations in the system for 
all energy densities so that all Class II operators not re¬ 
lated to energy conservation satisfy ETH, we obtain an 
analog of Eq. [TT] for the energy: E/Va > e max , and, 
(e ma xb — E)/Va > e max . Expressed in terms of the 
fraction Va/V , and the energy density of the eigenstate 
e = E/V, this constraint is 


/ < /* = min 



^max ^max 


( 20 ) 


Let us emphasize that the above constraint is a neces¬ 
sary condition for ETH to hold for all Class II operators, 
not a sufficient one. Just as some Class I operators can¬ 
not satisfy ETH when / is finite, we expect that there 
also exist Class II operators for which ETH fails when 
/ is finite, even when the above condition holds. Even 
so, significant deviation in the eigenvalue spectrum be¬ 
gins where this constraint breaks down, as our numerical 
results will demonstrate in Sec. ivm 
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Figure 1: Entanglement entropies Si through S 4 for a model with no conser vatio n law (left panel, given by Eq.pOlat L = 21), 
and a model with particle number conservation (right panel, given by Eq. 22 at L — 27 with filling N — 6). We use the 


parameters mentioned in the text to place each model at a nonintegrable point. In each case we consider eigenstates in the 
k = 1 sector, with subsystem size La = 4. The grey vertical line denotes infinite temperature (point of maximum Si), and the 
black circles mark the theoretical predictions for the entanglement entropies there. The brown markers denote the theoretical 
values of the entropies in the limit La, L 00 while La/L — y 0, as given by Eqs. 26 and[27| Notice that the Renyi entropies 


all match at infinite temperature if and only if there are no additional conservation laws besides energy. 


D. Summary 


III. A WARMUP: EIGENSTATES AT INFINITE 

T 


Let us summarize the discussion in this section. 


A. Von Neumann and Renyi entropy 


1. We conjecture that ETH holds for all local 
and non-local Class I operators as long as Va <C 

V. This implies that ETH is not restricted only to few- 
body operators (as can be seen in the limit Va/V -a 0 
as Va, V 00 ). When the subsystem is taken to be a 
finite fraction / < | of the total system size, we provide 
some evidence in Sec. |VII| that all operators not involving 
energy conservation satisfy ETH as well. 


2. We conjecture that ETH also holds for all 
Class II operators when Va V. It follows that a 
single eigenstate contains information about all energy 
densities available to the system. When / is finite, we 
suspect that ETH also holds for all operators that probe 
the system below an energy density e* (given by Eq. 19) 
and that do not involve energy conservation; however, we 
leave this as an open question. 


3. Determining the full Hamiltonian from a sin¬ 
gle eigenstate is equivalent to the satisfaction of 
Eq. [l]for both Class I and Class II operators. Our 

results provide strong evidence that this is true when 
Va <C V. Therefore, one should be able to extract infor¬ 
mation about the full Hamiltonian at arbitrary energy 
densities/temperatures using a single eigenstate. 


By definition, the thermal entropy reaches a maxi¬ 
mum at infinite temperature. Together with Eq. [35j 
this implies that when ETH holds, eigenstates at “in¬ 
finite temperature” are ones where the entanglement 
entropy is at its maximum. Consider a ID trans¬ 
verse field Ising model with longitudinal field, H = 
^2f =1 (ofof+i + h x af + h z cr?y Here the von Neumann 
entropy Si takes its maximum possible value when the 
eigenvalues of the reduced density matrix are all equiva¬ 
lent to one another. Thus, from counting the basis size 
of the reduced Hilbert space, we expect for infinite tem¬ 
perature eigenstates that each eigenvalue of the reduced 
density matrix will approach 2~ La in the thermodynamic 
limit when / = La/L < i From this, it follows that the 
Renyi entropies at infinite temperature satisfy 

S a = L A log 2 , ( 21 ) 

that is, they are independent of Renyi index a. The 
left panel of Fig. [l] shows how the entropies Si through 
64 together match this predicted value at the infinite 
temperature point for a L = 21 system with periodic 
boundary conditions and subsystem size La = 4. In 
general, as L -a 00 the T = 00 entropy density is given 
by S a /L A = log 2 . 

Now let us instead consider a model with an additional 
conservation law, namely particle number conservation. 













Consider a ID chain of hardcore bosons 


and the Renyi entropies are given by 


h = -J 2 (^+1+f'&.V 2 +h.c.) 

* ; ( 22 ) 

+ 2_ J (Vnin i+ 1 + V niUi+2) 

i 

where ni = b\bi. We focus on this system with periodic 
boundary conditions at the non-integrable point t = V = 
1 and t' = V' = 0.96. This model was previously studied 
and shown to exhibit ETH in Refs. 1601611 

Due to particle number conservation, the reduced den¬ 
sity matrix from any pure state is block diagonal, with 
each block corresponding to some filling number Na of 
the subsystem A. The block of the reduced density ma¬ 
trix corresponding to filling Na is a d^ A x d^ A ma¬ 

trix, where djsr A = (^). At infinite temperature and for 
La/L < 7 j, the eigenvalues of pA must be equal to one an¬ 
other within a given block, but the eigenvalues in differ¬ 
ent blocks will be different: they are in fact proportional 
to (nZ% a a ) , the number of microstates consistent with 
such a configuration in subsystem A. Taking into account 
that tr (pa) = 1, one finds that each of the d^ A = (^) 

eigenvalues of p^ A) are given by X Na = (^I % a a )/(n)- 
The spectrum of pA we find for a single eigenstate (as 
shown in Fig. [2]) is in agreement with that of the thermal 
reduced density matrix pA,th(P = 0) studied in Ref. 01] 
consistent with ETH. 

With this, the von Neumann entropy at infinite tem¬ 
perature becomes 

Si = - ^ dN A ^N A log Xn a (23) 

n a 


= log (E d N A ^N A J > (24) 

where the sums over Na are restricted to subsystem par¬ 
ticle fillings Na that satisfy the constraint in Eq. [9j The 
above expressions are valid when La/L < 1 . 

Because the eigenvalues are non-uniform, the Renyi 
entropies S a at infinite temperature depend on the Renyi 
index <a, in contrast to an energy-only conserving model. 
The right panel of Fig. [l] shows how the actual values 
of Si through S 4 match those predicted by the above 
counting argument. 

For comparison, we also calculate S a analytically in 
the thermodynamic limit. For simplicity, we consider 
the limits, L,N,La 00 such that n = N/L is held 
constant, while La/L —>> 0. In these limits, one can eval¬ 
uate the expressions in Eq.[24]using Stirling’s approxima¬ 
tion log (ad) ~ x\og(x) — x. One finds that in the limits 
considered, S a receives contribution only from Na given 
by 


TV* = 


La 


i+a-i r 


(25) 


Thus, S a probes the system at filling N\/La = 
r, which is different than the actual filling n, 


unless a = 1 (which corresponds to the von Neumann 
entanglement entropy). This also immediately leads to 
expressions for Renyi and von Neumann entanglement 
entropies in the thermodynamic limit: 


a; 
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eigenvalue index 


S a /L A = -F log [n a + (1 - n) a ] (26) 

a — i 

and 


Si/La = ~ [nlog(n) + (1 - n) log(l - n)]. (27) 

We plot these values in Fig.[l]for comparison. Remark- 
ably, even with the small system sizes we can access, the 
difference between the exact finite size result (obtained 
by counting over all sectors) and the result valid in the 
thermodynamic limit is quite small. 

In the above derivation, it is also possible to relax the 
restriction La/L -a 0 as L -A oc. We then find that 
N\ is given by the solution to 


Figure 2: Eigenvalue spectrum of the reduced density ma¬ 
trix of an infinite temperature eigenstate, Pa(\L) p=o) f° r the 
hardcore boson model Eq. [22] with L = 27, La = 4, and fill¬ 
ing N = 6. The red lines plot the theoretical value of each 
eigenvalue in the thermodynamic limit, determined from the 
filling N A of the sector in which it lies. 


NZ = 


La 


1 + 


( 


n-fN* A /L A 



which reduces to Eq. [25] when f 0. 

Let us note a few things about this equation: 


( 28 ) 
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1 . When a = 1, the solution is = tiLa, regard¬ 
less of /. Thus, the von Neumann entropy always 
probes the system at its given filling, even when / 
is finite. Further analysis shows that Eq. [27] holds 
generally when f < \. 

2. When the system is at half filling (n = |), the 
solution is iVJ = ^La, regardless of / or a. 

3. When <a>l, 0</<^, and n ^ the filling 
fraction N\/La probed by the Renyi entropy S a 
actually depends on /. As a result, the Renyi en¬ 
tropies for a given La depend on /. This can be 
contrasted with the von Neumann entropy, which is 
independent of / as long as / < The right panel 
of Fig. [l] illustrates this nicely: the analytical f 0 
prediction for the von Neumann entropy (Eq. |27| ) 
matches the corresponding numerical result quite 
well, but the Renyi entropies differ significantly be¬ 
cause / = 4/27 is finite. 

We expect that analogous features hold true also for the 
model that conserves only energy, which we will discuss 
in the later sections. 


B. Subsystem energy variance 

Let us also consider the average subsystem filling vari¬ 
ance of the particle-number conserving system given by 
Eq.[22]at infinite temperature. While the average subsys¬ 
tem filling is given by (Na) = nL a — Nf, the variance in 
this quantity for a single eigenstate with / = La/ L < ~ 
is given by 

((N A -(N A )) 2 ) = L A (l-f)(l-n) J L-. (29) 

Although both the filling and its variance are propor¬ 
tional to La as expected, the variance includes an ad¬ 
ditional factor (1 — /), which causes it to be suppressed 
compared with the grand canonical ensemble when / is 
finite. In Sec. |vn] we will witness a similar suppres¬ 
sion of the subsystem energy variance when the condition 
La/L —»> 0 is relaxed. 


IV. MODEL HAMILTONIAN WITH ONLY 
ENERGY CONSERVATION 

To develop some understanding of the questions posed 
in the introduction, we study a finite ID quantum spin- 
1/2 chain with the following Hamiltonian: 

L 

H = Y (°f°f +1 + Ka* + h z (T Z i) (30) 

2 = 1 

We set h x = 0.9045 and h z = 0.8090 such that the model 
is far away from any integrable point, and is expected to 


satisfy ETH in the sense of Eq. [l] as shown in Ref. 

We use periodic boundary conditions throughout. 

We diagonalized the Hamiltonian in Eq. [30] for system 
sizes up to L = 21, obtaining all eigenvalues and eigen¬ 
states. As hinted earlier, to each eigenstate we assigned a 
temperature /3 -1 by finding the value f3 for which the en¬ 
ergy expectation value in the canonical ensemble matches 
the energy of the eigenstate: 

WM _ tr {He-P*) 


By definition, (3 = +oc for the ground state and /3 = — oc 
for the highest excited state. In practice, the range of 
available /? values on a finite size system is much smaller. 
With L = 21, for instance, the first excited state has 
/3 ~ 4.0, and the second-to-highest excited state has 
P ~ —0.6 (as determined from Eq. 31). It follows that 
eigenstates outside the range 4.0 > (3 > —0.6 will not 
appear fully thermal due to the large thermal correla¬ 
tion length expected at low temperatures. (This can be 
seen for instance in Fig. [3J where the finite size correc¬ 
tions to the linear scaling of the entanglement entropy 
become more prominent as temperature decreases.) An¬ 
other thing to consider is that the infinite temperature 
eigenstate |V ; )/ 3 =o completely random and contains no 
information about the Hamiltonian. In a finite size sys¬ 
tem, states near infinite temperature will also contain 
little information about the Hamiltonian and will there¬ 
fore be unable to predict properties of the system at other 
energy densities. As a result of these finite size consid¬ 
erations, we typically study values of f3 between 0.2 and 
0.5 in the remainder of this paper. 



Figure 3: Scaling of the von Neumann entanglement entropy 
Si with subsystem size for the L — 20 system given in Eq. [30] 
Up to [3 = 0.5, the scaling is linear for small La, suggesting 
that the states are volume-law and are thus likely to satisfy 
ETH. The /3 = 1.0 eigenstate, on the other hand, is clearly 
not linear, and is too close to the ground state at this system 
size to exhibit ETH. 













10 






Figure 4: The von Neumann entropy Si and Renyi entropies S 2 , S 3 , and S 4 for the system given in Eq. [30]with L — 21 and 
La — 4. Here, Za — tr a(g~^ H a ). The entropies of the reduced density matrix at each energy density agree remarkably with 
the the entropies calculated from the canonical ensemble, given by Eqs. |2a| and |2b| 


V. VON NEUMANN AND RENYI ENTROPY 
OF EIGENSTATES AT FINITE T 

A. ETH prediction for von Neumann and Renyi 
entropies 


has an effective temperature /3 -1 . Z( l,/3) is the regular 
partition function of the system 3 Therefore, keeping 
terms only to the leading order in the subsystem size , the 
above expression leads to Eq. [3] advertised in the Intro¬ 
duction, 


Let us consider the Renyi Entropy S a = 
— log(tr PaH^} p)) corresponding to an eigen¬ 
state |^)/3 at inverse temperature /3. Assuming that 
ETH, as encoded in Eq.[2a| holds, S a may be reexpressed 
as: 


S 


a 


1 


a — 1 


log 


( Z(A,a,l3) \ 

V z {\,py ) 


(32) 


where Z(A,a,/3 ) is the partition function of the system 
on an a-sheeted Riemann surface, such that subsystem A 
has an effective temperature (a/3) -1 while subsystem A 


1 / e - a pv A f( a p)- a pv x f(p)\ 

Sa = “a - 1 l0g \ e -a/3VA/(/3)-a/3V^/(/3) J ( 33 ) 

= —V A p U{aP)-m) (34) 

a — 1 

where / is the free energy density. Therefore, the wave- 
function at temperature /3 _1 can be used to calculate the 
free energy at temperature (aft -1 . Indeed, the same re¬ 
sult also follows using the approximate form in Eq. |2b| 
Taking the limit a 1 leads to the conclusion that the 
von Neumann entanglement entropy Si satisfies 


ft = V A sM 


(35) 
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Figure 5: Scaling of the entropy deviation A S a = 
S a (pA, - S a (pA(\^) /0 )) with 1/L for constant La av¬ 
eraged over all eigenstates in the range 0.28 < ft < 0.32, for 
Si (top panel) and S 2 (bottom panel). The error bars rep¬ 
resent one standard deviation away from the mean. For Si 
this deviation is strictly non-negative, but for higher Renyi 
entropies it can oscillate and become negative before tending 
to zero as L -A 00 . 


Figure 6: Scaling of the von Neumann entropy deviation ASi 
with 1/L for constant ratio La /L averaged over all eigenstates 
in the range 0.28 < /3 < 0.32. As in Fig. [5j the error bars 
represent one standard deviation away from the mean. Even 
though this plot considers the case where the subsystem size 
La becomes infinite as L -A 00 , the entropy deviations are 
going to zero rapidly as L becomes larger. 



1/L 


where s t h(/3) = Si(pA,th(P))/L a is the thermal entropy 
density at temperature f3~ 1 . 


Figure 7: Scaling of the Renyi entropy deviation A S 2 with 
1/L for constant ratio La/L averaged over all eigenstates in 
the range 0.28 < /3 < 0.32. As in Fig. [5] the error bars 
represent one standard deviation away from the mean. 


B. Numerical Results for von Neumann and Renyi 
entropies 


Fig.[3]shows the scaling of von Neumann entropy Si as 
a function of subsystem size La for the eigens tates 
of our model (Eq. |3Q[). As discussed in Sec. IIC 1[ we 


expect Eq. [35] to hold as long as Va < Vj, in the limit 
Va,V 4 -A 00. This implies that in the thermodynamic 
limit, the function Si(Va) is expected to form an inverted 
triangle shape, similar to the behavior of a random pure 


state (Eq. [6|. However, in a finite total system at any 
non-infinite temperature, Si is an analytic function of 
the ratio Va/V with a negative sign for as shown 

in Fig. [ 3 ] (note that the sign of the curvature is fixed 
by the strong subadditivity of entanglement). However, 
even in finite system, the volume law does hold to a good 
accuracy when Va < H/2, and the finite size scaling, dis¬ 
cussed below, indicates that the inverted triangle shape 
is recovered in the thermodynamic limit. 
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Fig. [4] shows the comparison of Si, S 2 , S 3 , and S 4 cal¬ 
culated for each individual eigenstate for a subsystem size 
La = 4inaL = 21 system, with their ETH predicted 
canonical counterparts, Eqs. [35] and [3] We use two dif¬ 
ferent canonical counterparts corresponding to Eqs. [2a] 
and [2b] the latter version being susceptible to bound¬ 
ary errors, which nevertheless are expected to vanish as 
Va, Vj 00 • The agreement for each entropy is remark¬ 
able. It is worth re-iterating that the Renyi entropies 
for an eigenstate \^)p encode the free energy densities at 
temperatures different than /3 _1 (Eq. [3]), and these re¬ 
sults provide an instance of non-local Class II operators 
satisfying ETH. Also note that as a becomes larger, finite 
size effects become more pronounced because S a probes 
the system at lower temperatures (<a/3) -1 . 

We also studied finite-size scaling of the von Neumann 
entropy and Renyi entropies by keeping La constant 
and varying the total system size. The top panel of 
Fig. [ 5 ] shows the deviation — SthC#) for 

eigenstates in a range of temperatures. The difference 
ASi/La seemingly goes to zero faster than any inverse 
power of L, and is consistent with an exponential depen¬ 


dence ASi/La 


a —L 


, or at the very least, a power-law 


decay ASi/La ~ 1 /L X with x 1 (although we should 
caution that inferring the precise asymptotic finite size 
scaling behavior using exact diagonalization studies is an 
inherently difficult task). The bottom panel shows a sim¬ 
ilar plot for the deviation of Renyi entropy S 2 from its 
ETH predicted value, Eq. [3] The finite size scaling of 
AS 2 is relatively difficult because unlike Si , S 2 shows os¬ 
cillations as a function of La (see e.g! 41 * 63 *). Despite this, 
AS 2 is less than a few percent of S 2 itself. 

Fig.[ 6 ]plots the entropy deviation ASi/La for constant 
ratio La/L at all available system sizes. Although it 
is difficult to do a detailed scaling analysis with so few 
points, the data strongly suggests that ASi/La vanishes 
in the thermodynamic limit. 

The finite size scaling of Renyi entropies at constant 
ratio La/L is less conclusive, as can be seen in Fig. [7] The 
analytical argument for the particle number conserving 
model suggests that the Renyi entropies S a for a 7 ^ 1 
do not match their canonical counterparts when Va/V is 
held fixed. Ref. [64] arrived at similar conclusions using a 
different approach. 


VI. EXTRACTING THE HAMILTONIAN FROM 
A SINGLE EIGENSTATE 


In this section we will present numerical results that 
substantiate our conjecture that ETH is valid for all Class 
I and Class II operators when Va V as V 00 . 
Our numerical results consider the case where Va is held 
constant as V —> 00 . We expect that all results in this 
section also hold true when the limits are take n such 
that / = Va/V —>> 0 as Va, V 00 . In Sec. VII we will 


explore more carefully the case when / < \ is finite. 

We begin by probing in detail the entanglement spec¬ 


tra of individual eigenstates as well as the correspond¬ 
ing Schmidt states. Specifically, we compare four dif¬ 
ferent quantities, as shown in Fig. [8] which test the va¬ 
lidity of Eqs. [2a] and [2b] The agreement of the spec- 
trum of =±\og[p A (\il))p)} with that of =± log[pA,th(/?)] 
as well as with the actual Hamiltonian Ha in region 
A implies that essentially, the Schmidt eigenvalues A i 
satisfy A^ oc e~ f3EA ’ i where Ea,% are the eigenvalues of 
Ha • Similarly, the agreement with the expectation value 
(ui\HA\ui) shows that the Schmidt eigenvectors |iq) have 
the same character as the eigenvectors of the thermal 
density matrix. 

To probe the Schmidt eigenvectors further, we directly 
calculated the overlaps between the eigenvectors of the 
reduced density matrix Pa{\^)p) and the eigenvectors of 
the thermal density matrix pA,th(P) (see Fig. [9|. Again, 
we find excellent agreement. 

To quantify the extent to which Eq. [2a] is valid, we 
calculate the trace norm distance \ \pa(\^)/ 3) — pA,th(/3)\ |i 



Figure 8 : Comparison of the four quantities defined in the 
inset for an La = 4 subsystem at L — 21 and (3 — 0.3. Each 
quantity has been normalized so that the y-axis has units of 
energy density. The blue markers show the spectrum of the 
canonical (i.e. thermal) reduced density matrix while the red 
diamond markers correspond to the eigenvalues of a reduced 
density matrix pa {\^) /3 ) for a single eigenstate at tempera¬ 
ture /3; the grey markers show the eigenvalues of Ha with 
a shift ca = ^ log Za — ^ log trA(e _/3ifA ) so that it can be 
directly compared with — \og[pA(\^) p)] i n accordance with 
Eq. [2b] (note also that the combination Ha+ca is independent 
of the shift of the spectrum of Ha by an arbitrary uniform 
constant). Finally, the orange markers represent the expec¬ 
tation value of Ha , again with a shift ca, with respect to 
the Schmidt eigenvector \ui) oi p a{\^) p)■ In each case, the 
eigenvalues/eigenvectors are ordered from smallest to largest 
energy density. The horizontal lines plot the energy density e 
(dashed, grey) and the critical energy density e* = jL (solid, 
brown) of the original eigenstate \^) p, with respect to the 
ground state energy density of Ha + ca (dotted, black). 
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Figure 9: Overlap between the Schmidt eigenvectors | Uj) and the eigenvectors | (fi) of the canonical density matrix, for an 
L — 21 system with /3 = 0.3, and subsystem sizes La — 2,3,4, 5. In each case, the eigenvectors are ordered from most 
significant (largest eigenvalue) to least significant (smallest eigenvalue). 


between the reduced and canonical density matrices at 
various system sizes. The trace norm distance, defined 
as 


\\Pa(H) p)~PA,th{P)\\l = 


\/(Pa(\tP)p) - Pa, th(P)f 


(36) 

places an upper bound on the probability difference that 
could result from any quantum measurement on the two 
density matrice^. As such, it provides an excellent mea¬ 
sure of how distinguishable the two density matrices are. 
If the trace norm distance between two finite sized den¬ 
sity matrices is zero, they are equal to each other element 
by element. 


If ETH holds for all operators in subsystem A, then 
the results of Ref. [15] imply that the trace norm distance 
should go to zero as 1/L. The suggestion that the trace 
norm distance between the pure state and thermal re¬ 
duced density matrices with fixed subsystem size would 


tend to zero was also made in Ref. 66l We restrict our¬ 
selves to states in a /3 range given by 0.28 < (3 < 0.32. In 
the left panel of Fig. [lOj we plot the trace norm distance 
of every eigenstate in this /? range at La = 5 for a few 
select system sizes. For each system size, the distribution 
of the trace norm distance is nearly constant throughout 
the given /3 range. The right panel then takes this data 
for each pair of L and La and plots the mean and stan¬ 
dard deviation of the trace norm distance against 1/L. 
The trace norm distance is tending toward zero at least 
linearly with 1/L, perhaps even faster. 

These results, taken together, strongly support the 
conjecture that ETH, as given by Eq. [l] holds for all op¬ 
erators when Va <C V. The Schmidt eigenvalues and 
eigenvectors match at all energy densities, not just the 
energy density of the eigenstate. Our results also imply 
that when Va <C V, ETH as specified by Eq. [2a] holds. 
One consequence of this is that if Va is held fixed, the 
density matrices Pa{\^)p) and pA,th(P) become element- 
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Figure 10: Trace norm distance between the canonical density matrix pA,th(/3) and the reduced density matrix pa(\J>) /0 ) for all 
eigenstates | ip) p in the range 0.28 < f3 < 0.32. The left panel plots the trace norm distance for all such eigenstates with system 
sizes L — 12, 15, and 20, and subsystem size La — 5. The right panel plots the mean and standard deviation of the trace norm 
distance in this J range for values of L up to 21 and La up to 5. 


wise equal in any basis as F oo. 


VII. ETH WITH FINITE RATIO V A /V 


In this section we will consider to what extent ETH 
is valid when the ratio f = Va/V < \ is held fixed and 
finite as Va, V oo. As demonstrated in Sec. |VB[ the 


von Neumann entropy of pa(\^P) matches the thermal 
entropy in the thermodynamic limit even for finite f < \. 
In the current section we consider the extent to which 
other quantities match between a single eigenstate and 
the canonical ensemble. 

There is one notable Class I operator for which ETH 
(in the sense of Eq. [l]) fails when / is finite. As mentioned 
in Sec. |II C 1[ the subsystem energy variance taken from 
a single eigenstate is suppressed by a factor of (1 — /) 
compared with its value in the canonical ensemble. To 
demonstrate this relationship (given by Eq. [7]) , the top 
panel of Fig. [Tl] shows scaling of the subsystem energy 
variance with subsystem size La for both a single eigen¬ 
state and the canonical ensemble. While the energy vari¬ 
ance grows linearly for both for small La in both cases, 
the single eigenstate energy variance has an additional 
term that is negative and quadratic in the subsystem 


size. 


The bottom panel of Fig. [Tl] shows, for comparison, 
the variance of a different operator Ja between a single 
eigenstate and the canonical ensemble. The operator Ja 
(defined in the figure’s caption) is chosen to span the 
length of subsystem A and to include the same terms 
as Ha; however, the coefficient of each term is different. 
The fact that the variance of Ja matches between a single 
eigenstate and the thermal ensemble strongly suggests 


that all Class I operators that do not explicitly involve 
energy conservation will satisfy ETH in the sense of Eq. 
[TJ even when Va/V is finite. 

Let us now consider an implication of the difference 
in subsystem energy variance between Pa(\^)b) and 
pA,th(P)- This difference, which occurs only when / is 
finite, suggests that the trace norm distance between the 
two density matrices vanishes only when f 0. To 
explore this more carefully, note that the trace norm dis¬ 
tance places a bound on the difference in expectation 
value of any operator A that is bounded between zero 
and one as 


|tr(/>A) — tr(<rA)| < \\p — <t||i. (37) 

In order to calculate a lower bound on trace norm dis¬ 
tance due to the variance difference, we must write the 
energy variance as a bounded operator that maximally 
differs between the two density matrices. Naively one 
might be tempted to consider the operator Oa,p/ A 2 = 
(Ha — (Ha)^) 2 /AV, where for the operator to be 
bounded, A must be chosen to be the largest energy avail¬ 
able to the system. Since both O a,p and A scale linearly 
with V, the expectation value of this operator is actually 
zero in the thermodynamic limit for both Pa(\^)p) and 
pA,th(/3). Thus, no bound can be placed on the trace 
norm distance due to this particular operator. 

Let us instead now consider a modified energy variance 
operator, 


Aa,/3,a = Pa, p, a Pa, p, a, (38) 

where A is an arbitrary energy scale and Pa,p, a projects 
onto the subspace where Oa,p /A 2 has eigenvalues in the 
range [0,1], thus making A A,p,A a bounded operator. 
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Figure 11: (Top panel) Subsystem energy variance with re¬ 
spect to subsystem size La for both the canonical ensemble 
(blue circular markers) and a single eigenstate 1^)^ (red dia¬ 
mond markers), at L — 18 and (3 — 0.3. The inset shows the 
ratio between the energy variances at each subsystem size, 
which is expected to fit 1 — La/L in the thermodynamic 
limit (Eq. [7]). (Bottom panel) The variance of an operator 

Ja = J2t=.+ h^a-) + J2f= i _1 which in- 

eludes the same terms as Ha but does not relate to energy 
conservation, is plotted for comparison. Here, the quantities 
hx\ hx \ and are each taken from the uniform distribu¬ 
tion [—1,1]. 


This operator considers the energy variance within a re¬ 
stricted window of width 2A about the mean energy. 

To arrive at an approximate bound due to this opera¬ 
tor, let us assume that the energy histograms of pA,th{/3) 
and Pa{\^)p) are given by normal distributions with vari¬ 
ance cr^ h and = (1 — /)cr^ h , respectively. Since both 
distributions have the same mean, the difference in ex- 


Figure 12: Trace norm distance between the canonical density 
matrix and reduced density matrix for constant ratio La/L 
and 0.28 < /3 < 0.32. As in Fig. [lOj the error bars represent 
one standard deviation away from the mean. The horizon¬ 
tal lines indicate the approximate theoretical minimum each 
trace norm distance can take, based on the suppressed energy 
variance given by maximizing Eq. |39| 


pectation values is expected to be 


D = tr[pA,th((3)A-A,0,A] - tr[p y 4(|V’) /3 )A A , / 3,z 


i r A 
W-A 


o~E 2 /2a] 


e ~ E2 / 2 


^thV^7T 0~/p \J‘2n 


E 2 dE. (39) 


Given cr t h and /, it is possible to find A numerically such 
that D is maximized. Although A is proportional to y/V , 
the value of D itself is independent of V as F -A oo, 
since cr t h also scales with y/V. The maximum quantity 
D then provides a lower bound on the trace norm dis¬ 
tance between pA,th(P) and pa(\^P)p) in the thermody¬ 
namic limit 67 . 

Let us now turn to our results on the scaling of trace 
norm distance with system size when the ratio / = La/L 
is held fixed as L, La oo, which are shown in Fig. [l2j 
Although there are few points available for each ratio, the 
trend is clearly for the trace norm distance to decrease 
as L increases. The horizontal, dotted lines denote the 
theoretical minimum each trace norm distance can taken, 
given by Eq. [39j Remarkably, for each subsystem ratio, 
the trace norm distance rapidly approaches this lower 
bound, suggesting that the bound may actually provide 
the result in the thermodynamic limit. This in turn im¬ 
plies that other operators which do not involve energy 
conservation are likely to have equal expectation values 
for a single eigenstate and for the canonical ensemble. 

We now turn to results on the entanglement spectrum 
when / is a significant fraction of the total system size. 
As discussed in Sec. |II C 2[ if the constraint in Eq. [20] is 
violated, the entanglement spectrum cannot match above 
a critical energy density e* = e/f (see Eq. 19), where e 
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Figure 13: Comparison of the four quantities defined in the inset of Fig. [8] for eigenstates of an L = 21 system with La = 7 at 
/3 = 0.2, 0.3, 0.4, and 0.5. Each inset plots a 12-bin histogram of the log of the density of states versus the energy density: the 
solid blue curve from a single eigenstate Pa{\^) p) and the dotted cyan curve from the canonical ensemble pA,th(fi). We notice 
that in each of the four plots, the eigenvalues of the reduced density matrix corresponding to a single eigenstate (red diamond 
markers) begin to deviate significantly from the other markers (in particular, the eigenvalues of the thermal reduced density 
matrix i.e. the blue markers), as the energy density reaches the critical value e* (denoted by the solid brown line), indicating 
breakdown of ETH beyond e*. 


is the energy density of the state \^)p. Fig. 13 shows the 


comparison of spectra of four different quantities consid¬ 
ered in Sec.[VT]for several different energy densities of the 
reference state \pj)p with / = 1/3, at four different val¬ 


ues of f3. With / = 1/3, the energy constraint Eq. 20 


violated, and therefore we expect that the entanglement 
spectrum should deviate from the actual spectrum of the 
Hamiltonian at least beyond the critical energy density 
e* = e/f. We find for each value of ft that significant 
deviation starts to occur essentially right at this critical 
energy density. 


Surprisingly, even though the entanglement spectrum 
does not match the actual spectrum beyond the en¬ 
ergy density e*, the expectation values {upH a\ui)/La 
continue to match the energy eigenvalues of the actual 


Hamiltonian! To understand this phenomenon better, 
we analyze the different terms in Eq. [lT| As argued in 
Sec. IIC2, the only way (ui\HA\ui) can exceed the total 
energy E of the eigenstate is, if the ^boundary term, 


-^boundary — ^ ^ \j ^ \ ® \-^AA\^j ) ® lw)’ (^0) 


scales with the total system size. We find that this is in¬ 
deed the case, as shown in Fig. [l4| In agreement with the 
general considerations in Sec. |II C 2 , the Schmidt eigen¬ 
values deviate from their ETH predicted values beyond 
e* (Fig. 13) and become considerably smaller. 

To summarize the results of this section, we provided 
evidence that ETH holds for all Class I operators not 
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Figure 14: Decomposition of the energy density corresponding to an eigenstate amongst the three terms in Eq. ^]for /3 = 0.2 
(left panel) and — 0.5 (right panel) at L = 21 and La = 7. The dotted magenta line marks the ground state of H^- As in 
Fig. [8] the solid brown line denotes the critical energy density e* for subsystem A. 


related to energy conservation. For Class II operators, 
there is a critical energy density e* beyond which ETH 
definitely fails, even though, surprisingly, the eigenvec¬ 
tors still seem to be correct. (We leave the understanding 
of this result for future studies.) 


VIII. AN APPLICATION: EQUAL-TIME 
CORRELATORS AS A FUNCTION OF 
TEMPERATURE FROM A SINGLE 
EIGENSTATE 


In the previous sections we provided evidence that a 
single eigenstate encodes the full Hamiltonian. As an ap¬ 
plication of this result, we now calculate correlation func¬ 
tions at arbitrary temperatures using a single eigenstate 
| ip)g. The basic idea is similar to the relation between 
the Renyi entropies and the free energy densities (Eq. [3]). 

In particular, consider the correlation function, 


{0(x)0(y))p tTl 


trA (^(|t/>)ff)Q(a;)C>(y)) 
tr.4 (PaQ^p)) 


(41) 


where x,y are located in subsystem A , away from the 
boundary. Using Eqs. [2a| [2b| to leading order in the sub¬ 
system size, (O(x)O(y)}f 3 , n equals the expectation value 
of the operator 0(x)0(y ) at a temperature (n/3) -1 . 

Fig. [15] shows the expectation values of local operators 
within subsystem A as a function of /3, as predicted from 
a single eigenstate at inverse temperature /3q (indicated 
by a yellow dot on the red curve). We choose operators 
that are as far away from the subsystem boundary as 
possible so as to minimize the finite size corrections. Even 
though the agreement with the canonical ensemble is not 
perfect, the qualitative trends and the numerical values 
match incredibly well, given the modest total system sizes 
to which we are restricted. These predicted correlators 
also undoubtedly suffer from corrections expected due to 
the conical singularity at the boundary of A in Eq. [41] 


IX. SUMMARY AND DISCUSSION 

In this paper, we analyzed the structure of reduced 
density matrices corresponding to the eigenstates of 
generic, non-integrable quantum systems. We argued 
that given an eigenstate \^)g with energy density e and 
a corresponding temperature /3 _1 , the reduced density 
matrix for a subsystem A is given by 

Pa{ W)p) = Pa, th(P) 

where 


PA, th(/3) 


tr ( e~P H ) 


if the condition Va V is satisfied. This means that 
for a fixed eigenstate | one can always extract the 
properties of the corresponding Hamiltonian at arbitrary 
energy densities by taking Va/V -a 0 as the limits 
Va, V -A oo are taken. Remarkably, even when Va/V 
(< 1/2) is taken to be fixed and finite, one can still ac¬ 
cess many properties of the underlying Hamiltonian for a 
range of energy densities in the interval described in Eq. 

m 

We also introduced the notion of “equithermal” (Class 
I) and “non-equithermal” (Class II) operators. In a 
canonical ensemble at temperature T, the expectation 
value of Class I operators depends only on the properties 
of the underlying Hamiltonian at temperature T, while 
the same is not true for Class II operators. Our results 
strongly suggest that all Class I operators not involving 
energy conservation satisfy Eq. [l] as long as Va < V/2. 
For Class II operators not related to energy conservation, 
Eq.[l]is seemingly again satisfied as long as the constraint 
in Eq. |2Q| holds. 

We also provided analytical results for the Renyi and 
von Neumann entropies of infinite temperature eigen¬ 
states of a particle number conserving model. These 
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Figure 15: Equal time correlators for an L = 21 system plotted against inverse temperature ft. The blue dots denote the 
expectation value with respect to each eigenstate, the dashed cyan curve plots the expectation value in the canonical ensemble, 
and the red curve plots the expectation value predicted from a single eigenstate at /3o = 0.3 (yellow dot) by raising the La = 4 
density matrix to the power P//3q and rescaling it to have unit trace. 


results substantiate our numerical results for the energy- 
only conserving model. In particular, we find that the 
von Neumann entanglement entropy for a subsystem of 
size Va equals the thermal entropy for that subsystem 
as long as Va < H/2, and therefore follows the so called 
‘Page curved at all non-zero temperatures, thus gener¬ 
alizing the original work of Page and otherP^, and in 

agreement with the recent work on large central charge 
Qprjy49H5f] 

In this paper we only considered contiguous subsys¬ 
tems. It seems reasonable to conjecture that Eq. [2a] con¬ 
tinues to hold as long as the support of operator O can 
be chosen to lie in a subsystem which is not necessarily 
contiguous and whose volume satisfies Va V. This en¬ 
compasses the expectation values of local operators such 
as (O(x)O(O)), where 0{x) is localized at point x and \x\ 
can be greater than L/2 (where L is the linear dimension 
of the total system). 


Let us mention some of the practical implications of 
our results. Firstly, the fact that a single eigenstate en¬ 
codes properties of the full Hamiltonian could potentially 
be a useful numerical tool. For example, one could imag¬ 
ine targeting a finite energy density eigenstate of a Hamil¬ 
tonian H by variationally minimizing the energy of the 
Hamiltonian (H — E) 2 with respect to trial wavefunc- 
tions. The techniques in this paper would then allow one 
to access thermal properties of the Hamiltonian without 
directly calculating the partition function, which could 
be extremely helpful for Hamiltonians that suffer from 
the sign problem. 

Secondly, owing to the recent progress in single atom 
imaging techniques in cold atomic systemd^, one can 
now access non-local operators experimentall)®^^. 
This potentially allows one to check some of our predic¬ 
tions pertaining to the violation of ETH in cold atomic 
systems. For example, one can perform a quantum 
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quench on a low entanglement state which would at suf¬ 
ficiently long times lead to a thermal state in the same 
sense as Eq. [2a| This in principle allows one to determine 
the underlying Hamiltonian of a cold atomic system by 
performing tomography on a small subsystem to obtain 
the corresponding reduced density matrix, and then tak¬ 
ing its logarithm. 

We conclude by posing a few questions and future di¬ 
rections. 

In this paper we extracted equal-time correlators at 
different temperatures using a single eigenstate. It will 
be interesting to see if a similar method also works for 
unequal time correlators at arbitrary temperatures. The 
main difference is that this requires calculating expres¬ 
sions such as Eq. m at an imaginary exponent, and es¬ 
timating the effects due to the conical singularity in this 
case requires further study. 

As mentioned above, we expect that all our discus¬ 
sion carries over to time-evolved states as well since such 
states are expected to also have thermal behavior at long 
times in the same sense as a single finite energy density 
eigenstate. If so, does the time scale for thermalization 
for a given operator (i.e. the time it takes for the expec¬ 
tation value of the operator to become equal to its canon¬ 
ical expectation value) depend on whether the operator 
is Class I (equithermal) or Class II (non-equithermal)? 

Another question concerns the subleading corrections 
to the entanglement entropy. One expects that there 


always exist subleading area-law contributions to the en¬ 
tanglement entropy (either von Neumann or Renyi) of 
a single eigenstate. Are these contributions also cap¬ 
tured correctly in the entanglement entropies calculated 
via a thermal reduced density matrix? Perhaps a more 
interesting question is whether the mutual information of 
two disjoint intervals (which cancels out both the volume 
law contribution and the area law contribution) takes the 
same value for a single eigenstate and its canonical coun¬ 
terpart. 
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